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Abstract 


The  exact  theory  of  incompressible  potential  flow  through  pump  im- 
pellers with  logarithmic  spiral  (constant  angle)  blades  is  extended  to  include 
radial  or  conical  turbomachines  with  small  variations  in  vane  angle  and 
passage  breadth.  Approximate  formulas  are  obtained  and  charts  given  for 
smooth  entry  and  performance  characteristics  of  pump  or  turbine  configura- 
tions when  the  solidity  is  somewhat  greater  than  unity  and  when  the  major 
blade  deviations  are  confined  tc  the  inlet  regions.  However,  the  results  of  the 
present  analysis  are  general  and  may  be  applied  to  blade  rows  with  reasonably 
arbitrary  camber  lines. 


I.  Introduction 

In  recent  years  there  has  been  some  attention  given  to  the  problem  of 
calculating  the  performance  of  compressor  or  pump  impellers  with  a finite 
number  of  blades.*’  “ By  and  large  these  efforts  consist  of  solving  field 
equations  either  by  numerical  methods  (i.e.  the  relaxation  technique)  or  by 
variations  of  stream  filament  theories.  If  compressibility  effects  must  be 
considered  it  appears  inevitable  that  such  methods  must  be  used,  particularly 
if  any  of  the  flow  details,  e.g.,  velocity  distributions,  are  desired.  Further- 
more, since  the  impeller  geometries  commonly  employed  are  rather  complex 
the  question  of  analysis  is  made  increasingly  difficult.  Although  it  is  possible 
to  obtain  velocity  and  pressure  fields  for  rather  arbitrary  configurations  by 
the  methods  proposed  in  Refs.  1 and  i,  the  labor  and  time  required  for 
solutions  are  prohibitive  for  most  engineering  situations.  Moreover,  these 
calculations  are  for  specific  cases  and  do  not  permit  generalizations  to  be 
made  from  them. 

It  was  precisely  because  of  these  limitations  that  the  earlier  efforts  of 

3-4 

Spannhake  and  Sorensen  were  devoted  to  obtaining  analytic  solutions  for 
impeller  geometries  of  somewhat  simplified  shape.  Spannhake  obtained  che 
incompressible,  nonviscous  head  flow-rate  characteristic  for  impellers  with 
straight  radial  blades  of  arbitrary  radius  ratio  with  constant  impeller  breadth. 
Later,  Sorensen  extended  the  theory  to  include  the  effect  of  constant  blade  or 
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stagger  angle,  bur.  was  unable  to  obtain  exact  solutions  in  closed  form  lor  the 
overall  performance  parameters.  A.  Busemann^  in  1928  finally  solved  this 
problem  and  obtained  expressions  for  the  head  flow-rate  equation  as  well  as 
for  the  "shockless  entry"  condition.  His  analysis  was  restricted  to  the 
potential  flow  in  "two-dimensional"  impellers  with  a finite  number  of  logarith- 
mic spiral  blades.  These  results  are  exact  insofar  as  the  results  are  given 
in  series  which  may  be  calculated  to  any  accuracy  desired.  Busemann* s 
theory  gives  a good  approximation  for  many  CttBeS  provided  the  flow  is  pre- 
dominantly radial  or  two-dimensional.  However,  most  impellers  are  not 
two-dimensional  nor  are  the  impeller  blades  strictly  logarithmic  spirals. 

It  is  to  this  last  problem  that  this  present  work  is  mainly  devoted.  Un- 
fortunately, a truly  exact  general  analytic  solution  would  involve  great 
difficulties,  so  that  approximations  must  be  ltlhde.  The  program  to  be  fol- 
lowed herein  consists  first  of  obtaining  a rttore  general  form  of  Busemann's 
results  by  a slightly  different  method.  Inasmuch  as  this  theory  is  applicable 
to  turbines  as  well  as  pumps  both  sets  of  coefficients  will  be  given.  It  will 
then  be  seen  that  approximations  ir.  the  spirit  of  the  thin  airfoil  theory  can  be 
made.  The  approximate  problem  is  then  formulated  and  the  results  are  given 
for  . • all  variations  in  blade  shape  and  impeller  breadth.  It  will  also  be 
shown  that  under  certain  conditions  these  results  apply  to  ftowa  over 
cones  so  that  the  present  results  may  form  the  starting  point  for  a theory  of 
mixed  flew  impellers. 

In  the  course  of  this  work,  approximate  formulas  may  be  obtained  easily 
for  the  case  of  high  solidity.  These  have  been  found  to  be  extremely  useful 
for  a rapid  estimate  of  the  overall  performance. 


II.  Formulation  of  the  Problem 


in  this  and  the  following  three  sections  we  work  out  the  potential  flow 
through  a circular  array  of  equally  spaced  logarithmic  spiral  blades.  The 
breadth  of  the  impeller  is  assumed  to  be  constant.  Under  these  conditions 
the  flow  is  two-dimensional  and  hence  the  methods  of  conformal  mapping  may 
be  used.  An  important  simplification  in  this  problem  is  that  the  flow  due  to 
the  rotation  of  the  blades  may  be  treated  separately  from  the  through  flow. 


In  this  sec- 
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i.eM  the  flow  responsible  for  net  discharge  through  the  runner, 
tion  the  properties  of  the  mapping  function  which  enable  these  problems  to  be 
solved  are  worked  out.  In  the  succeeding  sections  the  flow  arising  from  the 
rotation  of  the  blades  in  still  water  (displacement  flow)  and  the  through  flow- 
are  obtained. 


Mapping  Function 


. 5 


The  mapping  function  is 


z = 


’ w - w~ 

IN 

~Vt  - l/w 

o 

o 

w - 1 

w - 1 

L o j 

L o J 

2iY 


or. 


dz 
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(1) 


(2) 


It  is  easy  to  verify  that  Eq.(l)is  the  desired  function  since  when  w is  on  the 
unit  circle  Eq. (2)  gives  Arg  dz/z  = y which  is  precisely  the  equation  of  a 
logarithmic  spiral.  Moreover,  Eq.(l)is  a single  valued  function  if  |w|i  1 
and  the  point  wq  is  not  enclosed.  If  w is  enclosed,  once  the  argument  of 
z increases  by  2tt/N,  thus  moving  a point  in  one  blade  passage  to  a congruent 
pcint  in  the  next  adjacent  one.  Fig.  la.  The  point  wq  is  the  complex  ( 
constant  and  represents  the  origin  of  the  physical  plane.  The  point  w = oo 
corresponds  to  z = oo  also.  This  functiev.  i»  slightly  different  than  the  one 
usually  used  for  annular  or  cascade  flows. ^ The  end  points  of  the  blade  are 
given  by  the  roots  of  dz/dw  = 0.  If  the  point  corresponding  to  the  exit  tip  of 
the  blade  is  arbitrarily  chosen  as  unity,  then 


and 


w. 


1 


w , 


ei(2r  + 25  + u) 


i 6 

w = a e 
o 


a = sin  y / sin(  y + 6 ) 


/ ->  \ 


(4) 


where  corresponds  to  the  inlet  edge  of  the  blade.  From  the  above 

the  exit  radius  of  the  blade  t ^ = [ z^j  = 1,  thus  the  radius  ratio  of  the  im- 
peller is  given  as  r^/r?  = jzj  or  using  Eq.(i)one  obtains 


1 + cos  2‘ 


r 1 _ sin^2  Y + 8 ) ^ 

r ^ ~ sin  6 


2(  n - 6-Y 


sin  2 Y 


where 


N In  T^/r i 
"cosy 


In  case  a « 1,  the  relation 
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For  straight  radial  blades  y = 0 and  Eq. (5)  reduces  to 


ri  r - - 1 1 
7-  = .mj 


F rom  Eq.  4 


tan  Y =■ 


a sin  6 
l -a  cos  6 


from  which  Fig.  lb  may  be  constructed.  As  an  aid  in  determining  the 
constants  a and  6 the  following  expressions  are  useful: 


a2-l 


2 cos^  Y 

,shK-tsinr(Tf-i.-jc?1 . i 

L 2 cos  Y J 


a2-l 


= 4 co9Y  e 


o~-4  sin  y (it  - 6 1 
1'  cos  y 


may  be  used  as  an  approximation.  Equation(5b)is  plotted  in  Fig.  2. 
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III.  The  Displacement  Flow 


i 

1 


1 

3 


In  order  that  there  be  no  flow  through  the  surface  of  the  blade,  the 
velocity  component  of  the  flow  normal  to  the  blade  must  satisfy  the  condition 

- n • <c  x r 

where  n is  the  outward  normal  of  the  blade  surface.  Thus  the  normal 
vector  velocity  must  be 

V — co  cosy  ' + 

nz  b 

Let  F be  the  complex  potential  of  this  flow.  Then 


dF  _ dF  dz 
cTw  Hz  * Hfw  * 


and  dF /dz  is  the  conjugate  of  the  vector  velocity  V in  the  z-plane.  Since 
polar  coordinates  must  be  used  in  the  w- plane,  we  have  as  the  boundary  con- 
dition on  the  unit  circle 


V 

rw 


m cos  y 


-1(tt/2+y) 

C • 


(6a) 


With  the  aid  of  Eq.  (2)  one  obtains  finally 


rw 
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where 
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(6c) 


The  problem  is  now  to  find  the  corresponding  tangential  velocity 
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* ! C"  V 

•,  rL 


distribution  on  the  unit  circle  subject  to  the  additional  restriction  that  this 
velocity  vanish  at  infinity.  The  solution  may  be  obtained  in  any  number  of 
ways.  The  most  direct  consists  of  writing  the  complex  potential  F in  the 


form 


TT/,„i  - 
■*•  \ ••  / 


oo 

V- ' a -n 

/L  Anw 

n=o 


and  solving  for  the  coefficients  by  means  of  Eq.  6b  and  hence  obtaining  the 
conjugate  velocity  V^.  The  procedure  that  will  be  followed  herein  is 
essentially  the  same  as  that  outlined  above  and  differs  only  in  some  details. 

In  order  to  effect  a solution,  the  fact  will  be  used  that  if  the  normal 
component  of  the  velocity  about  the  unit  circle  can  be  written  as 


Vr  = f(e-1®)  +7(elS) 


then  the  tangential  velocity  is 


..  1 I"  -r  , iOx  ,,  -id.  I 

v0  = T |_T^e  ) ‘ f(e  ) J 


provided  that  the  velocity  should  vanish  at  infinity.  This  statement  is  essen- 
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tially  a variation  of  Milne -Thompson' s circle  theorem.  Thus  in  order  to 
obtain  the  conjugate  velocity  Vqw  all  that  must  be  done  is  to  arrange  Eq.(6b) 
in  ascending  and  descending  powers  of  e*®  and  then  to  apply  Eq.  (7).  In 
Eq.  (7),  VQ  is  frequently  spoken  of  as  the  "conjugate  function"  or  simply 
the  conjugate  of  V • 


Conjugation  of  the  Displacement  Flow 


Since  |wq| -j  i , the  brackets  in  Eq.  (6b)  may  be  expanded  by  the 

binomial  theorem.  If  the  terms  are  arranged  in  positive  and  negative  powers 
iO 

of  e*  one  obtains  for  the  boundary  condition  in  the  w-plane 
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where 


A = V’  r(m+k-q)r(m-q) 

^ — 1 fm+k)l  ml  |w  | *'m 
m = o'  ° 


The  following  terms  must  now  be  conjugated: 
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The  terms  and  5^  are  easily  conjugated  by  a straightforward  application 
of  Eq.  (7).  However,  it  is  evident  that  if  the  denominator  of  is  expanded 

by  the  binomial  theorem  then  S,  will  contain  terms  of  both  positive  and  nega- 

iO 

tive  powers  of  e . Thus  in  order  to  accomplish  the  conjugation  it  will  be 
necessary  to  separate  these  terms,  can  be  written  as 


. ViY 

>3  = — nr- 
1 w 
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oo 


n=  ] 


,i(k-n)  0 


— n 

w 
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+ conj.  term 


with  the  aid  of  the  geometrical  series.  Now  if  in  the  first  sum 


k-n  = r = C,  1,2,,#  for  n ^ k 


k-n  = - r = 1,2,3 


for  n > k 


then  Sj  may  be  written  in  four  separate  terms  as 
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is  now  expressed  in  four  sums  which,  taken  two  at  a time,  are  arranged 
in  ascending  and  descending  powers  of  e*®  as  required  by  Eq,  ( 7).  It  should 
be  noted  that  the  boundary  condition  on  the  blade  must  result  in  a flow  which 
has  no  net  source  strength  within  the  unit  circle  in  the  w-plane.  Consequently, 
the  terms  corresponding  to  the  index  r = 0 in  are  either  zero  or  have  a 
sum  of  zero.  Hence,  in  treating  this  index  if,  dropped. 

The  conjugation  rule  of  Eq.  (7)  is  now  applied  to  terms  Sj,  and  S.^ 
to  obtain 
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where  A^  is  given  by  Eq.  (9).  The  resultant  sum  of  these  terms  may  be  con- 
siderably simplified  by  the  use  of  the  addition  theorem  for  binomial  coefficients 


g.  -mmmeiarjriwKt*  y 


The  result  of  these  operations  is 
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where  Fj,  are  hyper  geometric  functions  defined  by 


F . = F (-  q,  - q;  1;  —-—5) » F,  = lim  [ F (-  q,  - q;  £ ; - '--O 
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and  F ic  giver,  by  the  series 


F (a,  b;  c;  z)  = 


r,(c) 


r(a)r(b) 


P(a+n)r 
1 nl  P (c  +n) 


>+n)  n 


(10b) 


Unfortunately,  there  seems  to  be  no  way  to  reduce  the  double  series  to  a 
single  series.  However,  this  equation  does  provide  the  displacement  flow 
velocity  distribution  around  the  impelTer  blades  and  in  any  specific  case  pre- 
sumably these  terms  could  be  evaluated.  Equation  (iOa)  represents  only  part 
of  the  complete  displacement  flow,  for  in  order  10  establish  lift  or  torque,  the 
Kutta  condition  must  be  prescribed  at  the  »xit  tip  for  a pump  or  at  the  inlet 
edge  (Wj)  for  a turbine.  Consequently,  in  determining  the  head-flow  rate 

characteristics  the  values  of  V_  and  V0  will  be  needed,  i.e.,  those 

¥lw  ¥2w 


1 n _ 

A ^ 


velocities  corresponding  to  the  exit  and  inlet  blade  tips.  Fortunately,  the 
factor  of  the  double  series  (Eqs.  (10))  vanishe  s itt  these  points.  L_/C t V ^ ^ 
and  Vj  denote  the  tangential  displacement  velocities  at  the  exit  and  inlet 
tips.  They  become 
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(12) 


Equations  (11)  are  Busemann’s  principle  results. 


IV.  The  Through  Flow 

In  this  section  the  flow  from  a source -vortex  through  the  stationary 
(nonrotating)  array  of  blades  is  worked  out.  Inasmuch  as  the  through  flow  for 
normal  pumping  and  turbining  is  slightly  different  they  will  be  discussed 
separate  ly. 

a.  Pump.  The  point  wq  represents  the  origin  in  the  z-plane,  and  all 
of  the  through  flow  is  supposed  to  originate  there.  In  the  w-plane  the  solu- 
tion of  this  flow  is  obtained  by  imaging  a source  vortex  at  wq  in  the  'unit 
circle.  This  flow  potential  is 

G(w)  = log (w  - wj  + .£.l0g(i-wj  (13) 

where  A - Q + i P , C and  P being  the  source  and  vortex  strength  respec- 
tively. In  order  to  see  the  significance  of  A let  w--*wq  so  that  z — > 0 , 
Then  from  Eq.  (1) 

. -N 

w - w — » const  x z 
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Hence,  near  the  origin  in  the  z-pla.ne 


A M TSJ 

G(z)  — *•  -r-  log  z‘  = (Q  iP  ) (log  |z|  + i Arg  z) 
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and  the  stream  function  is 
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The  angle  between  a radial  line  and  the  tangent  to  the  streamline  is  given  by 


tan  a , = r 


= - P/Q  , 


I ’IfSr  ' 


where  a is  positive  as  shown  in  Fig.  3.  Hence,  it  can  be  seen  that  the 
streamlines  issue  from  the  origin  on  logarithmic  spirals  inclined  at  an  angle 
a,  to  any  radius. 

The  velocities  on  the  unit  circle  may  be  found  by  dinerentiation  of  Eq. 
(13)  and  they  are  found  to  be 
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The  radial  velocity  is  zero  as  it  should  be  and  upon  reference  to  Fig.  lb 
it  can  be  seen  that  for  w = 1 


Q cos  Y 


(tan  Y * tan  a^) 


(15a) 


"j*.  a 

< > F 

f 

i 


and  similarly, 


Q cos  y . . . 

V = — (tan  y - tan  a ) . 
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(15b) 


b.  Turbine.  In  the  case  of  a pump  the  flow  issues  from  the  origin  and 
approaches  the  blades  with  a prewhirl  characterized  by  the  angle  a^.  How- 
ever, in  normal  turbine  operation  the  flow  progresses  radially  inwards  and 
approaches  the  impeller  with  a given  angle  positive  in  the  same  sense 
as  aj.  The  potential  for  this  flow  is  obtained  by  merely  replacing  Q by 
- Q in  Eq.  (13).  Thus, 


1 


•w,.4 


\ I 


tan  = P/Q 


for  a turbine  and  also 


Q cos  Y 


itan  Y - lan 


\ (tan  Y * tan  q2^  « 

tj  tt  £ 


(17a; 


(17b) 


for  the  velocities  at  the  exit  and  inlet  blade  lips. 

In  the  development  of  the  characteristic  head  equat.on  Iot  a turbine  we 
shall  need  to  use  an  expression  for  the  circulation  which  must  be  put  around 
the  blades  in  order  to  satisfy  the  Kutta  condition  at  the  inner  blade  edge 
(w  = Wj),  However,  this  flow  differs  from  that  of  a pump  inasmuch  as  it  must 
not  disturb  the  flow  conditions  at  infinity.  Consequently,  the  potential  for  this 
flow  in  the  w-plane  must  consist  of  a vortex  at  the  point  wq  with  an  image  at 
the  inverse  point  within  the  circle.  Thus  if  H is  this  potential, 

rT  i rT 

H(w)  ~ i log  (w-wQ) log(w  - l/w'0)  • \ 18) 

The  velocities  at  the  outer  and  inner  blade  tips  are  then  respectively 


T,  = 


rT  p2 - 2 cos  y 


(19a) 


p j + 2 cos  y 


where  I is  the  strength  of  this  vortex. 


(19b) 


V,  The  Head  Flow-Rate  Equation 


a.  Pump.  From  elementary  considerations  the  head  of  a pump  can  be 


shown  to  be 


U,  V - U.  V 
H = 2 u2  1 U1 
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whe  re 


V 


u 


are  average  values  of  the  peripheral  component  of  the 


absolute  velocities  in  the  direction  of  rotation  at  the  exit  and  inlet.  The  head 


is  usually  expressed  in  terms  of  the  dimensionless  coefficient 


* = H/U-f/g  , 


and  can  be  shown  to  be  equal  to 


ryri 


, 2 
2ir  M 


where  P^,  are  the  values  of  the  circulation  at  the  outside  and  inside  of 


the  impeller.  If  p^  is  the  strength  of  the  circulation  about  each  of  the  N 


blades  then  this  coefficient  is 


n p 


* = - 


z 

r 


Zrt  r£to 


As  yet  there  is  presumed  to  be  no  circulation  around  each  of  the  blades.  In 
order  to  fix  the  value  of  P^  the  Kutta  condition  will  be  imposed  at  the  exit 
tip  of  the  blade,  i.e.,  it  will  be  assumed  that  the  value  of  the  velocity  must  be 
finite  there.  In  the  circle  plane  this  condition  is  that  the  velocity  at  w - 1 
must  be  zero.  This  condition  may  be  satisfied  without  altering  the  boundary 
conditions  by  placing  a vortex  at  the  origin  in  the  w-plane.of  such  strength 
that 

+ V.  - = 0 . 


d2+  vt2‘  “2T 


Thus, 


+ = - 


N 


to 


(V, 


+ V 


is  taken  as  unity.  The  first  term  is  independent  of  flow  rate  and  is 


incorporated  into  the  constant  . The  second  term  depends  only  on  flow 


rate  where  Q is  the  discharge  per  passage.  In  terms  of  the  dimensionless 
flow  rate  coefficient 


N Q 


U- 


2it  to 
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the  head  equation  assumes  the  form 

- \!'o  + <?  Cj_j  (tan  Y - tan  g.j  ) 

for  forward  curved  blades  (Y  > 0}  and 


(20a) 


^ - tP  C (tan  Y-  tan  a,) 


(20b) 


for  backward  curved  blades  (Y<  0).  The  shut-off  head  coefficient  \jf  is 
found  from  Eq.  (11a)  to  be 


4 co  s Y 

IT 


*o  = CH 


[a  coir  1 exp  [ - T - 6 )] 

[Fi  + N (4  • 0 F*\  ui) 


where 


_ 2 cos  Y 2a  sin  6 cot  Y 

c _ = 7 *— 

r2  1 + a - 2a  cos  o 


(22) 


This  coefficient  represents  the  effectiveness  of  the  blades  in  guiding  the  flow 
whereas  represents  the  effectiveness  in  bringing  up  the  average  circum- 

ferential velocity  to  the  impeller  tip  speed.  It  is  evident  from  Eq.  (21)  and 
Fig.  lb  that  0 ^ Cpj  - 1 . Values  of  the  coefficient  are  plotted  in  Ref,  8 

and  a more  detailed  plot  is  given  herein  in  Fig.  5. 

For  particular  values  of  r^/r^,  Y and  N,  a,  q,  and  <r  may  be 
determined  either  from  Eq.  (5b)  or  Fig.  2.  The  value  of  the  functions  F^ 
and  F^  may  then  be  found  through  the  use  of  the  series  expansion  Eq.  (10b) 
in  order  to  obtain  ^ . Values  of  this  coefficient  are  also  plotted  in  Ref.  8 
for  various  values  of  r^/ i^,  N and  y*  However,  these  graphs,  although 
quite  useful,  are  not  too  accurate  and  since  in  the  region  of  major  interest 
relatively  simple  expressions  may  be  developed  for  the  shut-off  head  coef- 
ficient, the  following  paragraphs  will  be  devoted  to  obtaining  more  convenient 
formulas  and  charts. 

Calculation  of  values  of  \|r0 


It  is  clear  from  Fig.  2 that  if  the  number  of  blades  is  small  and  the 


t 


- 1 S- 


radius  ratio  r^/r^  near  unity,  then  the  parameter  a L is  large  and  con- 
sequently the  expansion  of  the  hyper eeometric  functions  F,,  F,  will  contain 

“ a fc- 

only  a few  terms  so  that  these  cases  may  be  easily  worked  out.  However,  the 
majority  of  cases  for  pump  and  turbine  applications  turn  out  to  be  in  the  range 
of  the  parameter  1 £ tr/Ztt  £ 2.  An  inspection  of  Fig.  2 shows  that  the  quar. 
tity  jwQ|  = a is  about  unity  for  the  sc  conditions.  This  situation  is  for- 
tunate because  the  hypergeometric  functions  F^,  permit  easy  expansions 

to  be  made  when  the  argument  is  near  unity.  In  particular,  the  expansion 

9 

around  the  point  1 of  the  hypergeometric  function  is 


* F{c-a,  c-b;  1+c-a-b; 1-x) 


(23) 


provided  none  of  a,  b,  c are  negative  integers. 

Let  a 1 - a,  then  the  two  sums  F^,  become 


f.  = JZJi+siaL 

ln(  l+q)  P ( l+q) 


[l -SL3. 
L a+a 


. si iL-Mlzsl  a2 


A - .nr  a 

q+q  2(l-q-q)(q+q) 


] 


_ sin^asuMr^  iliiiaLClJLii)  A1+q*q  \ l + lltaillfa)  A + •••]  (24a) 

'(2+a+a)  L 2+a+1i 


rsimr(q+q)  P(2+q+q) 


F = Si\wl  , T i + ji  A+  • • • 1 - sin^qsin"q  ilLLtaLCiiiai  Aq+q 

N r (lrq)  P(l+q)  [ 1-q-q  J tr  sin  it  (q+q)  P(l+q+q) 


I 1 + — A + • * * 
L 1+q+q 


f 24b) 


With  the  aid  of  Fig.  2,  F,  and  F,  may  be  simply  computed  for  any  given 

16  2 

case.  Again  as  is  evident  from  Fig.  2,  A3*  a - 1 becomes  quite  small  and 

nearly  independent  of  Y for  about  <t/2tt  > 1.2.  Since  Eqs.  (24)  depend 
largely  on  1/a^  , it  cam  be  seen  that  both  they  and  the  coefficient  * will 
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rapidly  approach  the  limiting  value  given  by  ?.  = 1.  If  then,  only  iirst  order 


terms  in  A are  retained 


, 2 Y sin  2 Y « 


*o  = CH 
max 


4 cc  s Y 


(2  cosY ) 


L P ( 1+a) 


I r(H3t3)  _ r , 4 1 

lr(l+q)r(l+q)  L NJ 


1+q) • ( l+< 


4 cos 
s N 


4cos  Y h P(l+q)  P(l  + q)  P(l+q+q) 


(25a) 


CH  = 1 ' 


4 cos  Y 


or  more  approximately 


I - 


/ 2y  sin  2Y  s 
e*PV  - — - TJ — H 


(2  cosy) 


P(l+q+< 


4 cos  y r(i+q)  r(i+q) 


u*w 


The  gamma  function  factor  may  be  directly  calculated  from  its  product  defi- 
nition, however  a convenient  expansion  for  it  is 


.CQ+a+ai — = exp 

P(Hq)P(l+q) 


6.  5797  - 19.  233  — 


4 ? 1 
+ 8.6584  -4-^-  (8  zos  y-1}+  '** 
N*  J 


4 C 0 Y 

Although  this  scries  diverges  for  — ^ -■■■  >1  it  is  rapidly  convergent  for 
values  of  this  parameter  less  than  about  a quarter.  For  Y = rr/2,  q is  real 
and  any  of  the  tables  of  the  gamma  function  may  be  used. ^ From  Eq.  (25b) 
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lt  can  be  seen  that  ^ Q ~ 1 » and  the  equality  sign  only  holds  for  the  case  of 
N = oo. 


Figure  6 is  a plot  of  the  limiting  or  maximum  value  of  figured 
from  Eq.  (25b)  vs.  the  blade  angle  f for  various  numbers  of  blades.  If  the 
geometrical  parameter  <r/2ir  is  less  than  about  unity,  more  accurate  values 
of  the  shut-off  head  coefficient  can  be  computed  from  Eq.  (25a),  However 
the  proportions  of  most  designs  are  such  that  Fig.  6 gives  a good  estimate  of 
this  parameter. 


Condition  of  shockless  entry 

The  head  flow-rate  characteristic  (Eqs.  (20))  is  a straight  line  and  in 
potential  flow  any  point  on  that  line  could  be  achieved  in  operation.  However, 
in  the  flow  of  a real  fluid,  large  viscous  effects  may  be  expected  to  occur  for 
highly  unfavorable  pressure  gradients  in  the  inlet  portions  of  the  impeller. 
For  forward  curved  blades,  an  infinite  velocity  always  occurs  at  the  blade 
inlet  edge,  however  for  backward  curved  blades  there  is  one  flow  rate  for 
which  finite  velocities  result  at  the  inlet  and  this  condition  is  variously  known 
as  "smooth  entry"  or  "shockless  entry".  For  shockless  entry,  then,  the 
velocity  is  finite  everywhere  on  the  blade  and  presumably  real  fluid  effects 
should  be  minimized  at  or  near  this  operating  point,  other  things  being  equal. 
Hence,  the  present  calculations  should  be  expected  to  give  the  best  design 
predictions  at  or  near  smooth  entry.  ^ 


In  the  w-plane  \’(9j)  = 0 for  shockless  entry  and  hence 


0 


is  the  condition  that  must  be  satisfied.  The  value  of  the  circulation  P can 

z 

be  expressed  in  terms  of  the  exit  velocities  so  that 


(V  , - > + (Vf  - V.  ) 

dl  a2  l2 


= 0 


is  the  requirement.  From  Eqs.  (15) 


V \ ■ qs^-L  <tanT  - lan  ai>  <77  + 7? 


and  from  Eqs.  (11) 
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rjvf 

r - a 

r ”=] 

L"o^Tj 

Lw  - J 

Hence, 

NQ 


q — q 

i*  r W 1 * 


(tan  Y_  tan  q 


■ w p * W *j 

i>*"w4i  hFrNV 

O J Lw  - 1 J 


(28a) 


As  before,  most  interest  is  usually  in  the  case  of  <r/ 2-ir  >1  or  a - 1 <<■  1 . 
In  this  event  Eq.  (28)  may  be  evaluated  with  the  use  of  Eqs.  (5c)  and  (24). 
Then, 


4 cos^Y 


F1-NF2 


N sin  tr  q sin  it  < 


ir  sin  it  (q+q)  P(l+q+q) 


ii 3l  (a2-l) 


Eliiaial  .*Lii 
r(i+q)  r(i+q) 


or  for  4 cos  y/N  <<  1 


F,-NF,.(»t,VT)  r(u<l)r(: 

1 C 6rT  p(l+q+q) 


4 cos 


(4  cos^P^pfiX^iT)^) 


Equation  28  can  now  be  put  into  the  more  usable  form 

2 2 r 2 

(p  (tan  Y-  tan  t^)  = - (1  + ■ C--f~*-)  (~) 

*w  &N  2 


(28b) 


with  an  error  of  less  than  2 percent  if  (t/2tt  >1.2  or  so,  since  the  next  term 
gives  an  additional  factor  of 


, , 8 cos'  y tr  - 4y  sin  Y \ 

1 +— N -e»p(-  -r__) 


which  amounts  to  about  a percent  at  most  for  the  case  considered  here. 

The  significance  of  Eq.  (28a)  is  brought  out  by  noting  that  for  the  simple 
infinite  vare  theory 
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ri 

(tan  y-  tan  a ) = - U-) 
eoo  1 r 2, 
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in  which  the  subscript  oo  denotes  an  infinite  number  of  blades.  Thus  we  have 


Jl_  (i  +±L . 

'<Wx  3 Nc 


Since  ¥Q-  1 it  is  evident  that  the  flow  rate  for  smooth  entry  always  occurs 
at  a flow  rate  greater  than  predicted  by  the  infinite  vane  theory.  This  correc- 
tion is  important  since  from  Fig.  6 typical  values  of  \|rQ  may  be  0.  65  to 

0.  8.  Thus  if  corrections  are  not  allowed  for,  errors  of  30  percent  may  be 
made. 

Figure  7 is  a plot  of  the  head  characteristic  equation  (Eq.  (20))  for 
forward  curved  and  backward  curved  vanes. 

b.  Turbine.  The  same  considerations  are  used  in  developing  the 
characteristics  for  turbines  as  for  a pump.  In  fact  if  the  flow-rate  is  suf- 
ficiently high  for  a pump  with  backward  curved  blades,  a brief  inspection  of 
Fig.  7 will  show  that  the  bead  is  negative  and  hence  power  is  being  taken  out 
of  the  flow.  However,  turbines  normally  operate  with  the  flow  directed 
radially  inwards,  and  the  above  situation  is  termed  "reverse  turbine"  opera- 
tion. In  the  condition  of  normal  operation,  then,  the  flow  is  inward  and  hence 
the  Kutta  condition  must  be  applied  at  the  inner  blade  edges. 


The  head  drop  across  the  impeller  will  be  the  same  as  for  a pump, 
namely 


¥ 2 — 

Utz  “ 

The  Kutta  condition  at  the  inner  blade  edge  (w,)  requires  that 

V.  + V,  + V_  = 0 . 
fcl  dl  ^1 

With  the  aid  of  Eqs.  (lib),  (17b)  and  (19b)  one  obtains  in  a manner  similar  to 
that  of  Eq.  (20a) 
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The  subscript  "oT"  refers  to  the  turbine  "shut-off  head". 

It  is  of  interest  to  note  that  for  values  ci  y where 

<p  ^ *oT/CH  (tan  T-  tan  a2) 

the  head  of  the  fluid  increases  in  passing  through  the  turbine  an-d  thus  we  have 
the  analog  to  the  reverse  turbine,  namely,  the  inward  flow  pump. 

As  in  the  case  of  the  pump,  most  turbine  configurations  have  values  of 
the  parameter  cr/2ir  >1.2  which  means  that  ami  and  0.  In  this 

event  can  tje  closely  approximated  by 
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This  expression  is  seen  to  be  nearly  the  same  as  that  of  Eq.  (28a)  so  that  one 
has  approximately 
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The  behavior  of  the  turbine  shut-off  coefficient  is  thus  similar  to  that  of  the 
flow- rate  for  shockless  entry  of  a pump. 

Figure  8 is  a plot  of  the  characteristic  equation  (Eq.  (30))  for  a turbine 
for  various  values  of  a^.  For  normal  operation  < 0 and  the  sense  of 
rotation  is  opposite  to  that  of  a pump  with  backward  curved  blades.  The  con- 
dition of  no  leaving  whirl  usually  represents  a useful  approximation  of  the 
design  operating  point  for  most  configurations.  This  condition  is  easily  ob- 
tained from  the  inlet  velocity  triangle  to  be 

— = - tan  a,  . 
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£hockle  s s entry 

Reaction  turbines,  in  general,  have  blade  shapes  which  differ  substan- 
tially from  that  of  a logarithmic  spiral  so  that  smooth  entry  calculations  based 
on  the  present  theory  are  not  too  meaningful  except  possibly  for  a pump  of  con- 
ventional design  run  as  a normal  turbine.  In  this  case,  the  calculation  pro- 
ceeds in  the  same  manner  as  that  leading  to  Eq.  (28a),  i.  e.t  the  condition 
V(0=O)  = O must  be  satisfied.  When  this  occurs,  the  velocity  at  both  the  outer 
and  inner  blade  edges  is  finite  and  it  is  clear  that  this  condition  corresponds 
exactly  to  that  of  a pump  operating  at  shoekless  entry  except  that  the  direction 
of  flow  and  rotation  are  reversed.  Since  the  Kutta  condition  is  satisfied  at 
both  points,  the  flow  is  reversible  and  the  shockless  entry  computations  for  a 
pump  can  be  applied  to  turbine  operation. 

Torque  characteristics 

Turbine  performance  is  not  usually  presented  in  terms  of  the  dimension- 
less head  and  flow  coefficients  y and  ? unless  the  speed  is  maintained 
constant.  If  the  speed  is  not  constant,  torque  is  given  as  a function  of  angular 
speed  with  either  the  head  or  flow-rate  kept  constant.  In  the  present  work  it 
is  most  convenient  to  hold  the  discharge  constant.  The  torque  can  then  be 
expressed  by  means  of  the  dimensionless  coefficient  X 


**  PA2°r2Vr,  = 


(33) 


where  here  p signifies  density.  In  the  case  of  a real  fluid  the  efficiency  is 
given  by 

T>  . M 

‘t  ~ 

This  torque  coefficient  is  to  be  distinguished  from  that  of  a pump  operated  at 
constant  speed  which  is  given  by 

% = y (f 

and  where  the  efficiency  is  given 
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The  torque  coefficient  can  be  expressed  as 
' X = C„  (tan  Y - tan  a,)  - 

/I 


voT 


(34) 


with  the  use  of  Eq.  (30).  Thus  1L  is  a straight  line  function  of  l/'f.  The 

i . 

intercept  of  — * 0 represents  the  locked  rotor  torque,  and  the  value  of  l/tp 
<p 

for  "X  = 0 represents  the  turbine  runaway  speed.  These  results  are  shown 
in  Fig.  9. 


VI.  Extension  to  N o n 1 o g a r i t h m i c Spiral  Shapes 

The  theory  as  presented  thus  far  is  limited  to  flow  in  strictly  two- 
dimensional  radial  turbomachines.  The  simplified  shapes  demanded  by  this 
theory  are  seldom  realized  in  actual  designs.  There  is  no  doubt  that  the 
effects  of  blade  angle  and  passage  breadth  variations,  in  addition  to  meridional 
curvature  can  cause  important  deviations  from  the  two-dimensional  theory. 

The  flow  through  an  arbitrary  cascade  of  airfoils  is  a related  but  simpler 
problem  and  it  has  received  great  attention  in  the  past  twenty  years.  Un- 
fortunately, it  is  quite  difficult  to  get  an  exact  analytic  solution  even  for  a 
circular  arc  profile  although  special  cases  may  be  worked  out  by  Garrick's 
procedure.1^  Thi'-  difficulty  was  partially  circumvented  by  employing  approxi- 
mate methods  similar  to  that  of  the  thin  airfoil  theory.  The  work  of  Rannie^ 

is  among  the  most  useful  of  these  except  that  it  is  restricted  to  circular  arc 

l* 

camber  lines.  Integrals  for  other  shapes  have  been  given  by  Pistolesi 
but  not  evaluated. 

There  seems  to  have  been  relatively  little  analytical  work  done  for  an 

arbitrary  radial  array  of  pump  blades.  Calculation  schemes  using  Ackeret's 

method  of  the  continuous  distribution  of  singularities  have  been  worked  out 

15 

for  pumping  configurations  by  Betz  and  Fldgge-Lotz,  but  this  method  does 
not  permit  general  analytic  solutions  to  be  obtained. 

In  the  sections  that  fellow  a theory  for  nonlogarithmic  spiral  blades  will 
be  outlined.  Like  the  work  of  Rannie  and  Pistolesi  it  is  based  on  the  ideas  of 
thin  airfoil  theory.  As  the  svarting  point  for  this  method  the  solution  of  a 
more  general  boundary  condition  than  Eq.  (6)  will  be  obtained.  These  new 


1 


f 

1 


t ’ 


T -rr  v 


*4; 


solutions  arc  thcr.  interpreted  from  the  point  of  view  of  perturbations  on  a 
logarithmic  spiral  and  correction  terms  for  the  head  and  flow-rate  equation 
are  obtained.  The  resulting  blade  shapes  can  be  quite  general  and  no  restric- 
tion to  any  particular  variation  i*>  necessary.  These  solutions  ere  also  em - 
bodied  in  further  approximations  to  account  for  the  non-uniform  breadth  of 
most  impellers.  Finally,  a simple  correction  for  nonradial,  i.e.,  conical 
flows  is  given. 

a.  Solution  for  any  Power  Boundary  Condition 

Suppose  that  on  the  surface  of  the  blade  (Section  II)  the  normal  velocity 
is  to  be 

jvj  = kcocosY  rm+*  (>5} 

where  r,  u>,  y have  the  same  meaning  as  before  and  k,  m are  real  constants. 
The  vector  velocity  is  then 


„ . ..  I im  i(tt/2  + y) 

V = kucos  Y z,  z,  e ' ' , 

nz  b 1 b' 

and  following  the  procedure  leading  to  Eq.  (6b)  the  radial  velocity  in  the  w 
plane  becomes 
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and  since  jz^j  is  given  by  Eq.  (6c)r  j zh|  *“Tt*  may  be  found  by  repla' 
2 by  m+2  to  get 
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Thus  the  radial  velocity  in  the  circle  plane  can  be  written  as 
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where 


q'  = fT<  (1  +e2U  ) . 


or 


N<  = 


m + 2 


N . 


(37a) 

(37b) 


The  geometrical  constants  o£  the  mapping  are  ttril  determined  by  N,  r^/r  ^ 
and  y . Equation. (36)  is  now  formally  the  same  as  (6b)  except  that  q 
has  been  replaced  by.  q1,  The  tangential  velocity  corresponding  to  boundary 
condition  Eq.  (36)  is  ifre  same  as  Eq.  (11)  except  that  q and  quantities  de- 
rived from  it  are  replaced  by  q'  • 1 Ivj  s 
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is  F<-q'.  -q';  i ; 
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Equations  (38)  are  the  exact  solution  of  the  potential  boundary  value 
problem  given  by  Eq.  (35). 


b.  Application  to  the  Through  Flow 

Before  the  typesof  approximations  to  be  made  are  discussed  it  is  illus- 
trative to  show  the  applicability  of  Eqs.  (38)  by  using  them  to  obtain  the 
solution  for  the  through  flow  of  an  impeller  with  logarithmic  spiral  blades. 

If  in  the  z-plane  there  is  a source  flow  issuing  from  the  origin  at  an 
angle  o. ^ to  the  radius*  the  velocity  along  a streamline  can  be  shown  to  be 


NQ 
2 it  r 
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This  velocity  has  a component  -V  normal  to  the  blade,  and  in  order  to  make 
the  blade  a streamline,  a flow  must  be  added  which  has  a component  that 

will  cancel  che  contribution  from  the  source.  If  Y is  the  angle  of  the  blade 
to  a radius  then  this  component  is 

V = NQ  3in(r'ai)  _ 
n Zn  ? cos  a ^ 

By  inspection  of  £a.  35  it  is  seen  that  m = -2  and  that 

k = jig.  sin(Y~<,i) 

2ttco  cosycosn^ 


Upon  substitution  of  these  values  into  Eq.  (38)  one  obtains 


cos  ycos  Qj 


which  is  the  same  as  that  given  by  Eq.  (15a)  (except  in  slightly  different  form) 
for  the  tangential  velocity  at  the  point  w = 1. 

Boundary  condition  for  nonlogarithmic  spiral 

Let  us  first  consider  a log-spiral  blade  of  constant  angle  y . If  there 
is  a source  flow  at  the  origin  in  the  physical  plane  given  by  the  potential 

— (l-i  tanaj)  log  z , 


then  the  velocity  component  parallel  to  the  blade  is 

V = NS  C°S(r~°l) 

Zn  r cos  a ^ 


V-3 


9) 


This  flow  component  is  everywhere  parallel  to  the  blade  surfaces,  and  the 
required  solution  is  given  correctly  by  Eqs.  (15).  Suppose,  however,  the 
blade  angle  is  locally  slightly  different  from  y (Fig.  10).  The  condition  of 
no  flow  through  the  surface  is  then  that 

i 

V = 
n 


V tan  t 


(40) 


I 

where  V is  the  perturbation  velocity  normal  to  the  skeleton  line  A angle  y , 
V is  given  by  the  above  and  < is  the  angular  difference 


€ = Y - T , 141) 

r 

y being  the  angle  of  the  perturbed  blade.  It  will  now  be  assumed  as  is 
customary  in  thin  airfoil  theory,  that  boundary  condition  Eq.  (40)  may  be 
applied  on  the  logarithmic  spiral.  For  simplicity  it  will  also  be  assumed 
that  € « 1 so  that  tan  t = t . With  these  assumptions  the  boundary  condition 
on  the  logarithmic  spiral  may  be  expressed  as 


v'  = -NQ 

n 2-tr  r 


ces(r-aj) 


COS  O - 


If  the  angular  perturbation  e is  expressed  as  a power  series  in  the 
radius  the  required  solution  can  be  readily  obtained  by  means  of  Eqs.  (38). 

The  present  interest  in  this  work  is  for  pumping  configurations  and  the 
function  r-  (r  ) will  accordingly  be  specialized  to  suit  this  end.  For  this  pur- 
pose it  is  convenient  to  express  « in  terms  of  the  circumferential  variation 
y (Fig.  10)..  From  the  figure, 

o'  = 9 - y/r  , 

and  , 

• de 

tan  Y = r . 


Thus 


Now  if  c <<  1 


tan  Y = tan  Y £ 1 ■ cot  Y^p  + ^ cot  Y J 


» / e. 

tan  Y = tan  Y+  e/cos  Y » 


« = cos  y 


If -S5J 


The  variation  y is  useful  in  laying  out  impeller  blade  designs.  The 
true  length  of  the  blade  is  customarily  developed  by  graphical  procedures, 
and  typical  blade  variations  are  shown  in  Fig.  11  where  it  is  seen  that  most 
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of  the  departure  from  a logarithmic  - spiral  occurs  in  the  inlet  regions.  This 
fact  is  also  important  because  it  means  that  the  blade  shape  perturbations 
are  mainly  subjected  to  the  oncoming  inlet  flew.  Hence  for  this  reason  the 
perturbation  boundary  conditions  are  based  upon  the  component  of  the  inlet 
flow  parallel  to  the  angle  Y . This  situation  is  different  from  cascade  theory 
in  which  the  vector  mean  of  the  inlet  and  outlet  velocities  is  used. 

Perturbations  representative  of  pump  design  are 

y = X(l-r)3  (44) 

and 

^ = - 3 \ ( 1 - r)2  , 


so  that  according  to  Eq.  43 


c 


= X cos  y 


(y  - 3r  + 2r2) 


(45) 


The  coefficient  X is  determined  by  specifying  « at  the  inlet  (since  « = 0 
at  r = 1). 


The  equation  for  the  boundary  condition  now  reads 
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(46) 


The  constants  k and  m in  Eq.  (35)  can  now  be  identified  for  the  three  terms 
of  Eq.  (46) 
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With  the  aid  of  Eq.  (38),  each  of  the  contributions  of  Eq.  (47)  to  the  tangential 
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velocity  V may  be  determined.  They  are: 
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in  which  the  special  notation 
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has  been  introduced.  These  equations  may  be  evaluated  by  means  of  Eqs. 
(24)  and  in  the  case  of  <r/ Ztt  > 1.2  further  simplification  can  be  made. 

The  through  flow  solution  then  consists  of  the  above  correction  terms 
plus  the  main  contribution  from  the  log- spiral  solution.  In  like  manner* 
correction  terms  for  the  displacement  flow  will  be  obtained. 

c.  Application  to  the  Displacement  Flow 

The  boundary  condition  on  the  perturbed  blade  is 
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Vn=ro3cosY  =rcocos(Y+€). 

The  parameter  c is  again  supposed  to  be  small  compared  to  unity,  and  as  in 
the  preceding  section  this  boundary  condition  is  applied  on  the  logarithmic 
spiral  of  Y = const.  With  these  assumptions 

V = rocos  Y - « r « sin  y 

n ' 

is  the  condition  to  be  satisfied  on  the  biade  surface.  However,  the  first  term 
consists  of  the  displacement  flow  already  solved,  the  second  term  thus 
corresponds  to  the  perturbation  which  accounts  for  the  non-uniform  angle, 
hence 

Vn=  - r<o  c sin  y (49) 


is  the  additional  term  which  must  be  determined.  If  the  specialized  shape  for 
pump  blading  is  used  (Eq.(44))  one  obtains  finally: 
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The  constants  k and  N are 
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k^  = - 2 X siny  cos  Y ; m = 2,  N = 


The  additional  tangential  velocity  contributions  are  again  found  by  use  of 
Eas;  (38)  tc  be 
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The  effect  of  these  perturbations  on  the  head  flow-rate  equation  will 
now  be  considered. 

d.  Head  Flow-Rate  Equation 

From  Section  III  the  dimensionless  head  coefficient  is  seen  to  be 


where  the  primed  quantities  refer  to  the  correction  terms  given  by  Eqs. 

(48)  and  (52).  In  order  to  evaluate  these  expressions  the  condition  o-/2tt>  1.  2 
will  be  imposed.  Then  if 
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With  these  simolificatiens  the  head  eouation  becomes 
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Equations  (53)  and  (54)  represent  approximations  which  can  be  made 
if  the  solidity  is  sufficiently  high.  For  the  terms  ^(EN)  the  approximations 
are  satisfactory  but  are  less  accurate  for  f (N/E).  Although  more  accurate 
formulas  may  be  obtained  with  the  use  of  Eqs.  (E4),  the  present  results  pro- 
vide a fir  at -order  correc  ticn  to  the  head  flow  -rate  equation. 

These  corrections  are  applicable  for  pumping  configurations  in  which 
the  major  blade  snaps  variations  occur  in  the  inlet  portions  of  the  impeller. 

It  has  already  been  stated  that  the  shut-off  coefficient  and  the  guiding 
parameter  Cu  are  very  insensitive  functions  of  the  solidity  if  it  exceeds 
unity.  Consequently,  as  a result  of  inlet  perturbations  one  should  expect 
only  negligible  changes  in  performance.  Such  indeed  is  the  case  in  Eq.  (54) 
since  the  terms  are  all  of  the  same  order  of  magnitude  and  hence  almost 
cancel  each  other.  Thus  changes  in  inlet  angle  should  not  result  in  signifi- 
cant performance  changes.  However,  there  should  be  a change  in  the  shock- 
less entry  condition  at  least  of  the  same  order  of  magnitude  as  the  variation 
of  the  inlet  angle  itself. 

e.  Shockless  Entry 

The  additional  Contributions  of  Eqs.  (48)  and  (5E)  are  used  in  the  same 
way  as  the  development  leading  to  Eq.  (28a).  Thus  the  condition  for  shock- 
less entry  is 


(V.  - V.  ) + (V.  - V.  ) + (V.  ' - V.  ')  + (V, '-  V. ')  = 0 
tt  t2'  d,  d2'  tj  t,'  dj  d2' 


If  the  same  approximations  are  to  be  made  herein  as  in  Eq.  (54),  then  one 
may  use  the  result 
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The  shockless  entry  condition  then  takes  the  form 
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The  ratio  is  the  ratio  of  flow  rates  for  shockless  entry  with  and  with- 

out the  inlet  angle  perturbation.  A close  inspection  of  Eqs.  (55b,  c)  shows 
that 

«P^e  < 1 for  X > 0 

<Pc/<Pe  >1  for  X.  < 0 . 

This  behavior  is  to  be  expected  since  if  \ > 0,  Y > Y an<*  a "flatter"  blade 
inlet  angle  results.  This  can  be  easily  seen  if  the  limiting  case  N — » oo  is 
considered,  since  then 
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which  is  the  simple  infinite  vane  theory.  Equations  (55)  give  the  correct 
limiting  value  for  shockless  entry,  but  this  result  may  be  used  as  a "rule  of 
thumb"  for  determining  the  flow  rate  for  this  condition.  More  exact  values 
for  a finite  number  of  blades  can,  of  course,  be  computed. 
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VII.  Extension  to  Nonradial  Profiles 

Of  particular  interest  to  the  designer,  is  the  so-called  "mixed-flow" 
impeller  in  which  the  flow  approaches  axially  and  progresses  through  the 
impeller  on  more  or  less  conical  surfaces.  In  certain  cases  it  is  possible 
to  apply  the  results  of  the  foregoing  calculations  to  such  impeller  geometries 
and  thereby  obtain  information  in  this  important  field  of  application. 

The  principal  difficulty  in  such  mixed  flows  is  that  the  stream  function 
and  velocity  potential  both  do  not  satisfy  Laplace's  equation,  so  that  the 
methods  of  complex  variables  cannot  be  used.  This  circumstance  is  shown 
below  in  Eqs.  (56)  which  are  the  equations  to  be  satisfied  by  the  velocity 
potential  and  stream  function  T for  incompressible,  perfect  fluid  flow 
over  a conical  surface  with  varying  breadth  b and  constant  vortex  angle  £ 
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The  coordinate  system  is  shown  in  Fig.  12.  In  the  derivation  of  these  equa- 
tions it  has  been  assumed  that  the  curvature  of  the  wall  is  small  so  that 
variations  across  the  passage  breadth  may  be  neglected. 
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It  is  clear  from  Eq.  (56)  that  if  the  breadth  b is  not  constant  then 
and  Y both  do  not  satisfy  Laplace's  equation*  However,  an  important 
deduction  concerning  the  special  case  b = constant  can  be  made. 

a.  Flow  over  a Coni^cgil  Surface  of  Small  Constant  Breadth 

If  b is  constant  and  if  there  are  no  variations  across  the  breadth  then 
Eqs.  (56)  reduce  to  Laplace's  equation  in  the  conical  coordinates  , and 

the  two-dimensional  calculations  may  be  applied  as  follows:  First,  it  is 

clear  that  if  there  are  N blades,  that  the  flow  on  the  conical  surface  is 
periodic  with  a period  of  2ir  sin  K/N.  Thus  if  there  are  N blades  on  the 
cone,  they  are  equivalent  to  N/sin  K in  a straight  radial  array,  i.e. 

N*  = N/sin  K . (57a) 

The  other  geometric  properties  of  the  array  follow  in  a similar  fashion.  The 
radius  ratio  R^/Rj  = r^/r^  remains  unchanged.  The  solidity  parameter 
becomes 


$ . N1°S  r2/rl 

ZtT  T\ r sin  K cos  y 


(57b) 


The  boundary  conditions  for  the  displacement  and  through  flows  are  also  the 
same  as  for  the  pure  radial  tlow.  This  result  is  easily  shown,  for  on  the 
conical  surface 


Vn=  cii  x R cosy  = R « sin  £ cosy  = r u cosy  . 

Hence,  the  radial  flow  calculations  may  be  applied  to  the  flow  in  a conical 
impeller  of  constant  breadth  if  the  equivalent  number  of  radial  blades  is  given 
by  Eq.  (57a)  and  the  solidity  by  Eq.  (57b). 

An  interesting  limiting  case  is  provided  by  making  K — » 0 but  keeping 
crv  constant. 

If  is  the  slant  blade  height  then 
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is  a constant,  one  obtains 
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which  is  the  definition  of  solidity  for  a two-dimensional  cascade.  The  flow 
coefficients  of  the  radial  case  also  become  equal  to  the  cascade  results  since 
if  o'*  is  constant  and  N/sin  K — > oo  then  the  head  coefficient  (Eq.  (25a)) 
become  s 


which  is  the  cascade  result. 

b.  Effect  of  Variation  in  Breadth 

From  the  foregoing  discussion  it  is  evident  that  a realistic  treatment 
of  the  effect  of  breadth  variation  cannot  be  done  by  complex  variable  methods. 
However,  if  db/'dR  is  small  compared  to  unity  then  approximations  based  on 
complex  variable  methods  should  give  correction  terms  of  the  right  magnitude 
and  in  the  limiting  case  cf  an  infinite  number  of  vanes,  converge  to  the  exact 
answer. 

For  the  sake  of  simplicity,  in  this  section  it  will  be  assumed  that 
K - 90°,  i.e.,  a straight  radial  impeller.  In  Section  V-b,  the  through  flow 
solution  was  obtained  by  computing  the  "interference"  velocity  necessary  to 
make  the  blade  a streamline,  and  using  this  velocity  component  as  a boundary 
condition  in  Eqs.  (16).  Nov/  it  is  clear  that  if  the  breadth  is  varied,  the 
interference  velocity  will  be  changed  and  hence  a different  flow  will  result. 

In  this  section  the  additional  terms  for  the  through  flow  will  be  worked  out. 
Inasmuch  as  the  disp’^cement  flow  boundary  condition  is  independent  of  the 
breadth,  no  additional  correction  terms  should  be  expected  to  occur. 


The  boundary  condition 


V 
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NQ 

2irr 


sin(  v - a j ) b^ 
cos  a b(7) 


(57) 
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applied  along  the  blade  surface  makes  it  a streamline  in  the  through  flow 
where  b_  and  b(r)  are  the  impeller  widths  at  the  exit  and  at  any  radius 

i 

respectively,  and  y is  the  angle  of  the  perturbed  blade  as  before.  The  re^ 
striction  c = y*  ~ Y < < 1 will  be  retained  so  that 


sin(  y - q^)  = sin(  y - Qj)  + t ccs(y-flj) 


and  if  the  breadth  term  be  written  as 


b2/b(r)  = 1 + f(r) 
the  boundary  condition  becomes: 


NQsi^y-cij)  NQsin(y-a1) 


V = , 

n 2irr  cos  a 


smvY-Q.;  NQ  cos(  y - a . ) -i 

L f(r)  + * ±-  l + f(r) 

r cos  a | ' ' 2w  r cos  a. ^ L ' 'I 


The  first  and  third  term  of  this  expression  have  already  been  worked  out  in 
Sections  III  and  V-b  respectively.  Thus  the  additional  term  required  for 
breadth  variation  is 


_ NQ  sin(  Y-a,)  r n 

Vn  i,rcoia,  £<r)  [ 1 + « cot(T  ' “j)  ] • 


(59a) 


In  keeping  with  the  spirit  of  the  foregoing  approximations,  t will  be  neglected 
compared  to  unity.  It  should  be  pointed  out  that  this  assumption  results  in  a 
linear  theory.  In  order  to  obtain  the  solution  of  Eq.  (59)  by  means  of  Eq.(38) 
the  function  f(r)  must  be  expressed  in  a power  series.  Although  any  parti- 
cular function  could  be  worked  out,  it  seems  most  to  the  point  to  evaluate  the 
effect  of  breadth  variation  for  the  simplest  case  possible,  namely,  a linear 
one.  For  some  mixed  flow  designs  this  is  a fairly  realistic  choice  but 
Francis  type  impellers  require  more  complicated  functions.  Thus,  assume 


so  that 
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in  cos  a, 
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In  terms  of  £q.  (35)  one  has 


NQ  sin(  Y'  Cj)  , 

k.  = u. -» ; m=-2,  N - co 

1 r 2 it  03  cos  q j cos  y 


NQsin(Y-a.)  , 

k,  = - ja  -5 — — ; m - - 1 , N = 2N 

/ * f w t£\  r A a /«  rrt  2 v* 


:tt  co  cos  a ^ cos  y 


hence  the  perturbation  tangential  velocities  (found  from  Eq.  (35))  are 


Q sin(  Y- 


f F^NI-ZNF^ZN)  ,N 

p 2 2 cos  V 2' 


V = -p  i-  Jf 2N)  < — = + i--—rr  FJ2N) 
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Q sin(  y-  ttj) 
rr  cos  a . 


J(2N) 
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These  results  may  be  used  in  conjunction  with  Eqs.  (48)  and  (52)  to  obtain 
correction  terms  to  the  head  flow-rate  equation.  However,  one  may  expect  these 
terms  to  be  small  for  pumping  configurations  and  also  that  the  only  major 
change  will  be  in  the  condition  for  smooth  entry. 

With  the  former  restrictions  on  the  solidity  one  has  approximately 
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When  this  expression  is  included  in  the  shockless  entry  condition  Eq.  (55a) 

I 

remains  unchanged  except  that  A is  replaced  by  A , i.  e.  , 
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whe  re 


f r 


A = A + 


L ‘2 


(61) 


t» 

and  <?e  is  the  flow  rate  coefficient  for  shockless  entry  for  an  impeller  with 
variable  breadth  and  nonconstant  blade  angle.  The  coefficient  A is  given  by 
Eq,  (55b). 

If  p ; 0,  the  impeller  inlet  is  smaller  than  the  exit  (Eq.  (60)),  *o  that 
the  shockless  flow  rate  must  decrease  (Eq.  (60b)).  However,  the  usual  case 
is  that  the  iftlei  breadth  is  greater  than  that  of  the  exit  so  that  p.  < 0 and  the 
ratio 


tl 


In  the  limit  N — ► oo,  Eq.  (60b)  becomes 


tt 


(62) 


which  is  precisely  the  result  of  the  infinite  vane  theory.  If  there  is  no  pre 
whirl  (a^=  0)  this  limiting  condition  can  also  be  written  as 

tan  Yj  b2 
tan  r2  = ^ 


where  now  y Is  the  relative  flow  angle  at  the  inlet  for  two  different  breadths. 
Both  this  equation  and  Eq,  (62)  are  really  the  conditions  for  maintaining  the 
same  incidence  angle  to  the  blade  inlet.  Without  undertaking  detailed  calcula- 
tions this  result  is  probably  as  good  a rule  as  any  other  for  inlet  breadth 
variation. 

More  complicated  breadth  variations  could  be  used  to  account  for  dif- 
ferent profile  shapes.  However,  it  is  doubtful  whether  or  not  they  will  add 
significantly  to  the  results  of  Eq.  (60)  in  view  of  the  already  approximate 
nature  of  the  solution.  In  fact  because  of  the  obvious  crudeness  of  the  breadth 
approximations  we  have  ueen  somewhat  hesitant  to  present  them,  and  have 
done  so  only  because  we  are  unable  to  find  similar  computations  elsewhere. 
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c.  *^emar*  on  Application  to  Mixed  Flow  ^mpelle^rs 

With  the  aid  of  the  theory  developed  thus  far  it  is  possible  to  make  more 
extended  analytic  computations  of  mixed  flow  impeller  performance  than  here- 
tofore possible.  Thus  far,  the  effect  of  the  blade  forces  on  the  meridional 
streamline  distribution  have  been  neglected.  Tn  the  case  of  radial  or  mixed 
flow  pumps  the  deviations  can  probably  be  overlooked  at  least  at  the  design 
point  unless  blades  of  unusually  high  "warp"  are  employed  or  if  the  total 
head  distribution  is  not  uniform  across  the  passage  at  the  design  point.  In 
any  case,  these’effects  may  be  approximately  calculated  by  theories  recently 
made  available.*^' 


For  most  conventional  designs  with  constant  total  head  developed  across 
the  passage  it  will  still  be  necessary  to  obtain  the  meridional  streamline  and 
velocity  distribution  and  thus  work  out  a "pseudo -three -dimensional"  solution 
by  application  of  the  foregoing  theory  to  each  stream  lamina.  The  blade  sec- 
tions can  then  be  designed  to  have  smooth  entry  occur  simultaneously  across 
the  inlet  or  any  other  condition. 

This  discussion:  will  be  closedby  observing-  that  although  the  present 

theory  is  applicable  for  the  straight  radial  or  conical  profiles  it  is  still 

inadequate  for  the  smooth  entry  calculations  of  Francis  type  impellers  with 

rather  large  meridion  curvatures.  In  this  situation  the  flov  may  be  imagined 

to  take  place  on  an  infinitesimally  thin  lamina  of  revolution  and  the  methods 

IS 

of  conformal  mapping  can  be  apolied  to  this  "two-dimensional"  problem. 


d.  Remark  on  the  Influence  of  Compressibility 


A general  treatment  of  the  radial  flow  turbomachine  problem  is  too 
difficult  to  be  attempted  directly.  Specific  cases  have  been  worked  out  numeri- 
cally for  the  pressure  and  velocity  distributions  in  the  passage.  Radial  flow 
compressor  configurations  in  general  have  rather  high  values  of  the  solidity, 
e»g.,  <r/Zir  k Z are  not  uncommon.  For  these  c^ses  it  seems  reasonable  to 


suppose  that  the  design  coefficients 
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ivi  xoexua  uu  ut«  oi  two  or  three  since  the  relative  velocities 
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are  still  quite  subsonic.  That  this  is  indeed  the  case  is  shown  in  Ref.  1 
wherein  the  value  of  the  coefficient  ^ differs  by  less  than  one  percent  for 
an  impeller  operating  with  a tip  Mach  number  of  1,  5 and  zero.  In  fact,  the 
value  of  these  coefficients  for  compressible  flow  is'only  about  one  percent 
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dif£erent  from  the  results  given  in  Fig.  6. 

Thus,  it  would  appear  that  for  configurations  of  high  solidity  the  in- 
compressible calculations  give  a good  approximation  for  compressible  flow 
performance,  except  for  the  maximum  mass  flow  which  depends  mainly  upon 
the  inlet  design.  The  head  flow-rate  curve  should  then  be  nearly  the  same  if 
the  volume  flow  parameter  is  based  on  the  exit  volume,  i.e. , 

V,  . 

* _ Z u _ ^ 

* = ~ = p?7 

where  m is  the  mass  flow-rate/exit  area/sec  and  p,  is  the  exit  density. 
The  total  pressure  ratio  then  becomes 
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l 


CH  m (tan  y - tan  Oj) 


(63) 


where  k is  the  ratio  of  specific  heats,  J the  mechanical  heat  equivalent, 
PQ,  Tq  are  the  inlet  density  and  static  temperature  respectively. 


e.  Numerical  Examples 

In  this  section  a few  examples  illustrating  the  use  of  the  graphs  will  be 
worked  oiufc. 

1.  Consider  a radial  flow  pump  of  constant  breadth  with  6 blades  of 
constant  angle  y = 70°,  radius  ratio  0.  5,  and  with  no  pre-whirl,  i.e., 
a.j  = 0.  The  solidity  is 

'2—  - 6 In  2/2 it  cos  v = 1. 94 

so  that  Eqs,  (25b)  and  (29)  may  be  used  safely.  The  shut-off  coefficient  is 

(Fig.  6)  i|r  = 0.81  and  Cj.  = i.00.  Thus  the  head  flow-rate  relation 
T O H 

(Eq,  (2  0b))  for  backward  curved  blades  becomes 

* = 0.81  - 2.  745  0 , 
the  condition  for  shockless  entry  is  (Eq.  (29)) 


pmasf  , 


1 ,,  , 8tt^  cos^  70°x  _ , ,ft 

TT(i+,TT^ ) - 1*29 


oTST 


or  29%  more  than  the  simple  theory  would  indicate.  Now 


<Pe  = (r1/r2)2  cot  Y = 0.  091  , 


thus 


<p  = 0. 117 
e 


Since  there  is  r.o  inlet  whirl,  the  absolute  exit  flow  angle  is  given  by 


tan  a?  = */<p 


a2  = 76 


1 '5° 
1/2 


at  the  shockless  flow  rate. 


It  should  be  emphasized  that  ^ 0 is  only  a parameter  useful  in  deter- 
mining that  part  of  the  head  curve  near  tne  design  pvi  m,  urcsb  icil  fluid 
effects  prevent  ^ from  being  realized  near  shut-off.  • 

2.  If  it  is  desired  to  operate  the  pump  of  Ex.  1 at  a shocklesa  flow  rate 
of  <p  = 0.  100,  then  the  inlet  angle  Yj  will  have  to  be  increased  (i.e.  made 
flatter).  Thus  from  Eq.  (55a) 


^e  _ 0.  100 

077 


= 0.  855  , 


Since  y , N,  r^/r^  are  known,  X in  Eqs.  (55b,  c)  can  be  found.  For  g^=  0 

A = 0. 0324  X 
B = - 0. 224  X 


hence 


X = 0. 566  . 
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By  Eq  (45) 


t = 0.0664  = 3.  8° 


or  y^73.8°.  The  inlet  angle  must  be  increased  by  about  4°  in  order  to 
have  smooth  entry  occur  at  0 = 0.  100.  At  this  operating  condition  the  inci- 
dence angles  to  the  blade  is  5.  5°  whereas  in  the  case  of  Ex.  1 it  is  5°  at 
shockless.  For  this  case  then,  the  simple  rule  of  maintaining  a <~cnstant 
angle  of  attack  or  incidence  to  the  blade  will  satisfy  the  smooth  entry  con- 
dition. For  smaller  values  of  \/2tt,  however,  deviations  from  this  ruie  must 
be  expected  to  occur.  Without  further  calculation  we  may  also  suppose  tfc't 
the  same  result  holds  true  for  minor  breadth  variations. 


3.  Let  the  pump  of  Ex.  1 be  operated  now  as  a normal  turbine.  In 
general  there  will  be  some  whirl  leaving  the  impeller  so  that  a j / 0.  For  the 
normal  direction  of  operation  Eqs.  (30)  and  (32)  give 


= <p(2.  75  - tant^)  - 0.  322  . 

Typical  values  of  a.,  are  in  the  interval  -80°  £ a,  £ -60°.  If  the  value  of 

^ ^ 0 
is  chosen  is  that  for  the  normal  pump  operation  in  Ex.  7,  i.e,(  = - 78.  5 

then 

= 6,  9 1 9 -0.32 

is  the  performance  of  the  pump  in  Ex.  I operated  as  a turbine.  The  condition 
of  no  exit  whirl  is  the  same  as  that  of  the  shockless  flow  rate  in  Ex.  1,  viz. 

9 = 0.  117. 

4.  In  the  last  example  a 70  angle  pump  with  5 blades,  constructed  on  a 
45°  cone  with  a radius  ratio  of  0.  7 will  be  considered.  The  solidity  of  this 
configuration  is 


* N In  r,/r, 
or  _ 2 1 

2tt  ~ i-rr  sin£  co  s 


Ln  2 


Zn  sin45°  cos  70° 


= 1.  17 


which  is  large  enough  for  ihe  approximations  of  Eq.  (25b)  and  (29)  to  be  used. 

# . 

The  equivalent  number  of  radial  blades  is  N = N/sin  K = 7.  1,  hence  from 
F igs.  5 and  6 

CH  = 0.983 
= 0.  83 
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The  coefficients  i|f  given  in  Fig.  6 are  the  maximum  values  for  the  given 
number  of  blades.  Unless  C.^  falls  several  percent  less  than  unity  closer 
approximations  of  do  not  have  to  be  made.  In  that  event,  however, 

Eq.  (25a)  will  provide  a first  order  correction. 


For  0^=0 


, the  head  equation  is 
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If  the  breadth  of  the  impeller  is  comparable  to  the  radius  then  the 
analysis  will  have  to  be  done  for  a number  of  streamlines.  It  is  customary 
• to  design  pumps  in  such  a way  so  as  to  have  a constant  head  developed  across 

the  exit  breadth.  The  meridian  velocity  distribution  is  then  only  slightly 
. affected  by  the  blade  loading  and  one  may  use  any  means  to  obtain  the  meridian 

flow.  However,  at  "off  design"  flows,  the  effect  of  the  blade  loading  will 
become  important  for  small  values  of  K , i.e.»  shallow  conical  pumps,  and 
recalculation  of  the  meridian  velocity  profiles  by  the  methods  of  Ref.  16  may 
become  necessary. 


VIII.  Concluding,  Remarks 

There  have  been  three  purposes  in  presenting  this  material.  The  first 
has  been  to  outline  an  pvarf  theory  for  potential  flow  in  impellers  with  legs 
rithmic  spiral  vanes.  The  second  has  been  to  describe  an  analytic  "thin 
airfoil"  theory  which  could  account  for  the  effect  of  blade  angle  variation. 
This  approximate  theory  starts  from  the  solution  of  the  potential  flow  through 
a radial  array  of  log-spiral  blades  with  the  normal  velocity  component  being 
prescribed  as  any  power  function  of  the  radius.  A method  is  shown  whereby 
these  solutions  may  he  used  to  obtain  correction  terms  for  small  variation? 
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in  vane  angle  and  passage  breadth  and  also  when  the  meridional  streamlines 
are  cones. 

The  last  objective  has  been  to  develop  the  theory  in  a form  suitable  for 
radial  flew  pump  and  txirbine  configurations.  Formulas  and  charts,  useful 
for  rapid  performance  estimates,  are  provided  for  this  field  subject  to  the 
restriction  that  the  solidity  be  greater  than  about  1.2.  However,  general  re- 
sults are  also  given  so  that  other  designs  may  determined. 

It  would  have  been  desirable  to  include  a number  of  detailed  impeller 
computations  in  which  account  is  taken  of  the  meridional  streamline  distribu- 
tion. The  work  involved  is  considerable,  however,  and  for  that  reason  such 
results  are  not  presented  here. 

It  would  also  be  of  interest  to  extend  the  present  theory  to  account  for 
the  effect  of  strong  meridional  curvature  (i.e.  Francis  type  impellers).  This 
would  complete  the  theory  of  the  "pse'ido-three-dimensional"  impeller  with 
a finite  number  of  blades.  The  writer  is  interested  in  this  problem  and  would 
like  to  see  examples  of  such  calculations  made  available  so  that  the  im- 
portance of  these  parameters  could  be  determined. 
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List  of  Symbols 


IX.  Notation 


A - Area 

F - Complex  velocity  potential,  hyper  geo  me  trie  function  as  noted 
N - Number  of  impeller  blades 

Q - Flow  rate  through  one  impeller  passage 
T - Torque 

U - Circumferential  velocity  = ra 

V - Absolute  velocity 

W - Relative  velocity 

C„  - Flow-rate  correction  factor 
rl 

a - | wo| 

b - Breadth  of  impeller  passage 

g - Gravitational  constant 

k,  m,  n-  Constants 

r - Radius  in  physical  plane 

q - Complex  constant 

w - Complex  coordinate  in  circle  plane 

wq  - A constant  in  the  w-piane  representing  the  origin  of  the  z-plane  » ae 

z - Comnlex  coordinate  in  physical  plane  = x + ig 


q - Angle  between  absolute  flow  and  radius  vector, 
^Positive  increasing  counter-clockwise) 

X - Angle  between  blade  tangent  and  radius  vector 

5 - Argument  of  w^ 

« - Angular  perturbation  of  blade 

£ - Cone  angle 

- Angular  coordinates 

X - Constant 

a - Radiu?  vector  in  w-nlane 

I a 

Vln  r Jr 

tr  - Solidity  = 

cos  y 

- Torque  coefficient  at  constant  speed  - 

$ - Velocity  potential 

<p  - Flow-rate  coefficient  = NQ/ZitA^V^  -•  V /V? 

Od  - Torque  coefficient  at  constant  flow-rate  - $ 


Notation  (continued) 


- Stream  function 

- Head  coefficient  = H/U^/g 

- Angular  speed 


Superscripts 

- Conjugate  i.  e.  , z = x + i y,  z = x - iy 

- Conjugate  function 


Subscripts 

- refers  to  impeller  blade  quantities 

- displacement  flow  tangential  velocity  component 

- conditions  at  shockless  or  smooth  entry 

- normal  velocity  component 

- radial  "elocity  in  w-plane 

- through  flow  tangential  velocity  component 

- velocity  in  z-plane 

- circumferential  component  of  velocity  in  z-plane 

- .velocity  in  w-plane 


- tangential  velocity  in  w-plane 

- circulation  flow  velocity  component 

- conditions  at  no  flow  for  a pump 

- conditions  at  no  flow  for  a turbine 

- conditions  at  inlet 

- conditions  at  exit 

- conditions  based  on  an  infinite  number  of  vanes 
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Fig.  3 - Definition  sketch  showing 
direction  of  positive  a. 


Fig.  4 - Definition  sketch  showing 
inlet  anri  exit  veloc’tv  triangles  fo 
backward  curved  vanes. 
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Fig.  7 - Head-flow  rate  characteristic  line  for  a 
pump  impeller  operating  in  frictionless  incom- 
pressible fluid  with  a finite  number  of  blades. 


Fig.  8 - Theoretical  head-flow  rate 
characteristics  for  a turbine  rotor 
v/ith  a finite  number  of  blades  opera- 
ting at  constant  speed.  Note  that  the 
flow  coefficient  axis  is  the  negative 
of  that  for  a pump. 


Fig.  9 - Theoretical  torque  coef- 
ficient vs,  speed  coefficient  for  a 
turbine  operated  at  constant  flow 
rate. 
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Fig.  10  - Sketch  illustrating  para- 
meter for  a perturbed  blade  shape. 


Fig.  11  - Sketch  showing  type  of  per 
turbations  useful  in  pump  designs, 
i lie  plane  shown  is  the  logarithm  of 
the  physical  or  real  plane. 
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Fig.  12  - Definition  sketch  for  mixed  flow  conical  pump  of  variable  breadth 
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